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Many thermodynamic instabilities in one dimension (e.g. DNA thermal denaturation, wetting 
of interfaces) can be described in terms of simple models involving harmonic coupling between 
nearest neighbors and an asymmetric on-site potential with a repulsive core, a stable minimum 
and a flat top. The paper deals with the case of the Morse on-site potential, which can be treated 
exactly in the continuum limit. Analytical expressions for correlation functions are derived; they 
are shown to obey scaling; numerical transfer-integral values obtained for a discrete version of 
the model exhibit the same critical behavior. Furthermore, it is shown in detail that the onset 
, of the transition can be characterized by an entropic stabilization of an -otherwise unstable-, 

^ ■ nonlinear field configuration, a soliton-like domain wall (DW) with macroscopic energy content. 

O ' The statistical mechanics of the DW provides an exact estimate of the critical temperature for 

a wide range of the discretization parameter; this suggests that the transition can be accurately 
viewed as being "driven" by a nonlinear entity. 
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I. INTRODUCTION 



> 

In a seminal paper published a quarter-century ago, Krumhansl and Schrieffer [y explored the possibility that 
I ' nonlinear excitations (solitons) might drive structural phase transitions. They succeeded in identifying signatures of 
■ O "domain walls (DW)" in the thermodynamic properties of the one-dimensional 0^ model and gave a phenomenological 
^ [ account of a salient feature of structural phase transitions, the appearance of a central peak in the dynamical spectrum. 
However, important issues remained open. In particular, because the demonstration was limited to a system which 
does not exhibit a genuine (i.e. finite-temperature) thermodynamic phase transition, there was no possibility to 
explore the role of exact, or nearly exact, nonlinear field configurations in determining critical behavior, static or 
dynamic. 

A class of one-dimensional models, which has been proposed to describe a wide variety of thermodynamic insta- 
bilities, such as thermal DNA denaturation ^, wetting of interfaces and other similar phenomena, could 
[ prove to be of interest in the above context. Typically, such models include a short-range interaction between near 
. neighbors, and an on-site potential with a single minimum and a flat-top; in the context of DNA denaturation, as 
■ the temperature increases, particles begin to access the region of the flat top at an increasing rate, until, at a finite 
temperature, a macroscopic instability occurs, with a divergence in the average displacement (the "effective length" 
of the hydrogen bond which holds the double-stranded chain together), a divergent correlation length, and, generally, 
all the characteristics of a thermodynamic phase transition. In the simplest case, that of a harmonic interaction and a 
Morse on-site potential, the thermodynamics can be calculated exactly by using functional integral techniques |^ and 
S mapping the problem to the quantum mechanics of the Morse oscillator. The thermodynamic transition corresponds 

to the quantum mechanics of the disappearance of the last bound state. 
^ In this paper, we will show that soliton physics can offer a correct interpretation of such thermodynamic instabilities. 

O In particular, a relatively brief calculation can show us why and when a thermodynamic transition occurs. The logic 
can be summarized as follows: An exact, static nonlinear configuration of the continuum field can be found, which 
J> , "interpolates" between the low- and the high-temperature phase. We interpret this to a be a DW. The total energy 
of the configuration diverges, because each site which is in the high-temperature phase contributes a finite amount 
. ^ of bond (and elastic) energy. The infinity is in a sense a "minimal requirement": as long as the DW energy is finite, 
• no phase transition can occur in one dimension . A closer look at the properties of the DW is more revealing. In 
" " ' contrast to other Klein-Gordon field theories with a finite-energy kink, there is no zero-frequency Goldstone mode. It 
costs energy to shift the position of DW so that it can include more sites at the high-temperature phase. Conversely, 
there is an energy gain by "zipping back" to the low-temperature phase. At finite temperatures, the picture must be 
corrected to include entropic effects. The "phonon cloud" which accompanies the DW provides an entropic gain which 
ultimately balances the energy cost of extending the DW to macroscopic scale. The formation of a thermally stable 
DW is thus seen to "drive" the instability. At higher temperatures, it is entropically advantageous for a DW to extend 
itself by "unzipping" more sites towards the high temperature phase. The picture sketched above, which concentrates 
entirely on the thermal stability of the nonlinear interface, is more than a qualitative nonlinear scenario. It provides 
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an estimate of the transition temperature which (i) turns out to be exact within the continuum field framework and, 
(ii) can be nonperturbatively extended to describe lattice systems up to quite high levels of discretization. 

The paper is structured as follows: Section II provides a full account of exact results for the model. This is 
important to establish notation and context for the discussion. Moreover, although some results have appeared in 
the literature before, they are scattered in various sources; in addition, results for the correlation function and the 
susceptibility are new and the calculation has not been reported elsewhere [Q. Section III deals with the soliton 
physics of the transition. A brief discussion of results and further perspectives, with special emphasis to effects due 
to lattice discreteness, is given in Section IV. 



II. EXACT THERMODYNAMICS 
A. Model and notation 



The Hamiltonian of the model is 



Pn 

2m 



K 

+ iriVn- Vn-lY + y{yn) + Dhayn 



(2.1) 



where m is the reduced mass of a base pair, y„ denotes the stretching of the hydrogen bonds connecting the two bases 
of the n^^ pair and Pn = m{dyn/dt) corresponds to the conjugate momentum of y„. 

In addition to the kinetic and nearest-neighbor potential energy terms, ( |2.lD contains: (i) an on-site potential which 
describes the interaction of the two bases in a pair; the Morse potential 



Viyn) = D (e- 



(2.2) 



plotted in Fig. has been chosen because it has the correct qualitative shape; and (ii), a field-dependent term, which 
describes the effect of a transverse, external stress h. 

A set of parameters which has been used in the DNA denaturation context ^ is: a dissociation energy D = 0.03 
eV, a coupling constant K — 0.06 eV/A^, a spatial scale factor of the Morse potential a = 4.5 A^^, and a mass 
m — 300 a.m.u; the time scale is determined by = (K/my^^. The dimensionless ratio R — Da?' / K is traditionally 
used to distinguish between the "order-disorder" (discrete, i? ^ 1) and "displacive" (continuum, R<^1) regime; the 
given set of parameters corresponds to R = 10. L 

The classical thermodynamic properties of ( |2.l| ) can be described exactly in terms of the transfer integral (TI) 
equation (pi p| 



dy e 



-l3[V{x) + V(v)+2W{x,v)+Dah{y+x)]/2 



-I3e„ 



4>n{x) 



(2.3) 



where (3 = 1/(^^7"), fc^ is the Boltzmann constant and T is the temperature. 

In general, we are interested in situations for which h = 0; however, the external field is useful in practical 
calculations as a mathematical device. For example, by letting /i ^ 0+ it is possible to extract the scaling behavior near 
the transition; moreover, since the partition function is now divergence-free at any h > (and numerical calculations 
can in principle be performed at finite, decreasing h), previous criticism of the model on formal mathematical grounds 
lH is addressed. 

In the thermodynamic limit, the free energy per site (ex cept for a nonsingular term due to the integral over 
momenta) is determined by the smallest eigenvalue Eq of (2^), i.e.. 



f = £q In 

" 2/3 



27rm 



Other thermodynamic properties of interest are (i) the order parameter 

r'+oc 

a ^<y>-- 



dy y |0oP 

where (j)o is the normalized eigenstate corresponding to the eigenvalue Eq and (ii) the correlation function 



(2.4) 



(2.5) 



C{j) =< [y, -cj){yo-cj)> 



(2.6) 
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FIG. 1: The Mor se potential V{y) is shown (solid line) as a function of the dimensionless field ay. The groundstate wavefunction 
4'oi.y) (Eq. ( 2.19 )) is shown for different values of the temperature: dashed line for T — Tc/2.6, dotted line for T = Tc/l.S. 
Note the logarithmic scale for (jjoiu)- 



where the sum runs over states other than 0o, 

A„ EE /3(e„ - eo) , 
and the ofF-diagonal matrix elements M„ are given (in Dirac notation) by 

Mn =< n\{y- < y >)|0 >=< n\y\Q > 



(2.7) 
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B. The corresponding Morse oscillator problem 

In the gradient expansion (continuum) approximation, 

My) = M^) + <f>'nix){y + <(x)(y - x)V2 + 0{iy- xf) 



(2.9) 



which is strictly valid in the temperature window D <C ksT ^ D/R [Q, the integral equation (2.3) can be well 
approximated by the second-order differential equation 



1 



2f3^K dy^ 



Die-'^y -ly +Dhay 



4>n{y) 



(2.10) 



where e„ = e„ + ln(27r//3i4r)/2/3. The one-dimensional statistical mechanics TI problem is thus mapped to the 
zero-dimensional quantum mechanics problem of the Morse oscillator. 



In the following, we will list and/or derive some general properties of (2.12) at /i = 0. The change of variables 



z = 2dexp{-ay) with S = -V2DK , 

a 

the transformation (j)n{y) = e~^/^ z" Wn{z) and s = S^l — e„/L), leads to 

z „ -I- (25 -I- 1 - zj— hnw„ = , 



(2.11) 



(2.12) 



where n = S — s — 1/2. If n is a positive integer, the solution of ( 2.12| ) is a Laguerre polynomial 
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C. Bound state(s) and associated distinct length scales 



Noting that 4>n(y) remains finite over the interval [0, +oo[ only for positive values of s, we obtain the spectrum of 
bound states, 

^ = with n = 0,l,...,i?(<5-l/2), (2.13) 

where E{5 — 1/2) is the integer part of 5 — 1/2. It follows that as long as 5 exceeds a critical value 6c = 1/2, the 
ground state remains bounded. In the quantum problem, this provides a criterion for the critical mass of a particle 
below which it is driven out of the potential well by quantum fluctuations. We note that this is a general property of 
asymmet ric po tential wells; symmetric wells support a bound state for any value of 6 jl^, 19 1. 
Using ( ^.11 ), this Sc defines a critical temperature 

2V2KD 

Tc - —, 2.14 

in terms of which S = Tc/2T. As the temperature approaches Tc, the last bound state becomes less and less localized, 
and < y > in creases sharply, indicating the "melting" of the system. 
Using ( ^.4|) results in a free energy per site / = eo + /o, where 



and 

|2' 



eo = D 



I - t 



iiT<Tc (2.16) 



~ D otherwise 



where t = T/Tc — 1 is the reduced temperature; it follows that the entropy per site can be written as the sum of a 
non-singular part 

Snon-sing = S{Tc) + ks log (2-17) 

and a singular part, 

2D 

Ss^ng = if T<Ta (2.18) 

= otherwise. 



We note that (2. IS) implies a jump discontinuity of the specific heat at Tc; this is consistent with a specific heat 
critical exponent a ~ 0. 

For Tc/3 < T < Tc, the (normalized) ground state 



r(2(5- 1) 



g-./2 ^5-1/2 ^ (2.19) 



where F is the gamma function, is the only bound state. The asymptotic behavior (poiu) oc e "'"^ '^/^^v defines the 
spatial extent of the ground state 

The order parameter, obtained from Eqs. ( |2.5| ) and ( ^.19| ), is ||l^ 

cr = ^ [ln(2(5) - V(2(5- 1)] (2.21) 
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FIG. 2: Left scale: dependence of the order parameter a (diamonds) and its fluctuations (triangles) as a function of the 
redu ced t emperature t; the points were obtained numerically using the TI method; the solid line shows the analytica l esti mate. 
Eq. (2.22), with Sc = 0.36. Right scale: the susceptibility x as a function t (squares: numerical TI; solid line: Eq. ( ^.4(]| ) with 
Sc = 0.36). 



where ^p is the digamma function; in the vicinity of the critical temperature, this reduces to 



a{2S- 1) 2 2a5c 



(2.22) 



in the language of phase transitions, a oc \t\^ ^ where /3 = —1 is the critical exponent for the order parameter, i.e. the 
order parameter dz'^er^e. at the instabihty! Fig. | shows that the agreement of numerical TI values with ( p.22[ ) is 
excellent, provided one takes into account the numerical TI value of 5c — 0.36 (corresponding to Tc — 427K) for the 
discrete system. 

Local fluctuations of the order parameter are described by |l2| 



In the vicinity of the critical temperature, this reduces to 



A 



a(2(5-l) 2 



(2.23) 



(2.24) 



which implies oc 1^1^"^ with — —j3 — 1. Close to the transition, the order parameter and its fluctuations are 
therefore comparable, as shown by Fig. ^; in the wetting literature they are both interchangeably referred to as 
"transverse correlation lengths" . 

If we denote the lattice constant by I and identify the next-to-lowest eigenstate with the bottom of the continuum 
band, we obtain the inverse longitudinal correlation length 



^1 



(2.25) 



i.e. the corresponding critical exponent i/^ =2. The solid line plotted in Fig. (|^) shows a perfect agreement for the 
difference of the two first eigenvalues, between Eq. ( 2.25 ) and numerical TI results. 
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FIG. 3: Dep endence of (ei — eo)/D on the reduced temperature. The sohd Une corresponds to the exact (continuum) result 
cf. (2.25). The diamonds are numerical TI results. 



D. Static structure factor and susceptibility 

In terms of the correlation functions derived above, the static structure factor reads 

+ CXJ 



iqjl 



+ 00 +00 



j — — 00 K— 1 
+ OC 



sinh Ak 



cosh A„ — cos{q£) 



(2.26) 
(2.27) 
(2.28) 



where the sum extends over all the scattering states of Eq. ( 2.10| ). They are given by the second solution of 
the confluent hypergeometric equation ( 2.12 ) for imaginary s = in. In the notation of Eq. (13.1.1) of Ref. ||13|], the 
solution w{z) — U{l/2 ~ S + in,! + 2iK,z) has the correct asymptotic behavior: it grows as z^-^/^-"'-'^ as z ^ +oo 
(i.e., y —^ — oo), so that the exponential prefactor in the definition of (jj^, dominates the vanishing of the wavefunction 
at the Morse hard core. These scattering states correspond to plane wave superposition 



w{z) (xl-e2'^(''''^)z-2''' 



(2.29) 



as z — !■ (i.e., y +oo) where 0{k) ~ argr(l + 2iK) + argr(l/2 — d — in). The unnormalized continuum eigenstates 
are therefore given in the asymptotic limit y ^ oo by 



and the corresponding eigenvalues by 



D I 



As and e„ differ by a quantity that does not depend on k, from (2.31) and (2.7) we identify 



(2.30) 
(2.31) 

(2.32) 
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Near the critical instability and in the long- wavelength limit q£ ^ I, the static structure facto r can be calculated as 
follows for a given value of | . Due to the rapid oscillations of the scattering eigenstates ( ^.3(]| ) , matrix elements will 
vanish for large values of k whereas, for small k, the phase shift can be approximated by 6* = — arctan(KAa). Since 
t he g round state extends to large values of y, for most of the interval of integration we can use the asymptotic form 
(2.30) of the continuum states. Normalizing in a box of size {0{L^), L/2) to account for the soft core, the ground and 
excited states are respectively 



My) 



2i 



4>K{y) ^ sin(Kaj/ + 6'(k)), 



where L is the system size. We obtain 



/ 8A3 sin(cr + I 
M« ~ i\ — 



L 1 + K^A^a^ 

where a = 2arctan(KAa) = —29. This explicitly gives tan(0 + ct) = n\a and hence 

8A3 



(2.33) 
(2.34) 

(2.35) 
(2.36) 



In order to substitute the sum over k by an integ ral in Eq. (2.28), we need the correct density of space in k space; this 
can be obtained by se tting the eigenfunction of ( 2.34| ) equal to zero at the boundary x — L/2 and leads to a density 
of states aL/2'iT. Eq. (2.28) can be rewritten as 



^((Z,T) = A^^m||) 



(2.37) 



where 



Fig. 



F{x) 



cosh(arcsinh(a;)/2) 



demonstrates that the reduced quantity 

S{q,T) _ F{q£,\\] 



S{0,T) F(0) 



-4F(g^ll) 



(2.38) 



(2.39) 



obtained in the continuum approximation gives an excellent description of the discrete (TI) results over many decades 
of the dimensionless variable qS^; note that in the second equation we have made use of the property F{0) = 1/4. 
The zero stress isothermal susceptibility is given by 

h^Q \dh 
= f3Da^S{0,T) 

= —\t\-^ 



1 T 

25c Tc 



i.e. the corresponding critical exponent is 7 = 4. Again with 5c. = 0.36, Fig. 
(2.40) with TI values over many decades. 
At temperatures very near Tc and finite q, i.e. for q^|| S> 1, F{x) ^ 2/x'^ and hence 



(2.40) 

demonstrates excellent agreement of 



S{q,Tc) cx 



2-»7 



(2.41) 



with 77 = 0. 
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FIG. 4: The reduced static structur e fact or S{q, T) /S{0, T) as a function of for different values of the reduced temperature. The 
sohd hne is the theoretical curve (2.3!:) whereas the symbols were obtained using the TI method for different values of the 
temperature. 



E. The DNA order parameter 



The quantity which is directly accessible in DNA denaturation in the fraction of bound base pairs, which can be 
measured using UV absorbance. The probability P{y < b, T) of finding a given base pair at an equilibrium distance 
smaller than h (equal to the fraction of bound pairs, with a proper choice of 6), is given by 

Pi{y<b,T) ^ f dy \My)\^ (2.42) 

•/ — oo 

- 1 - (2.43) 

r(2(5-i) ^ ^ 

« {2S-1) Ei{25e-''^) (2.44) 

where 7 is the incomplete gamma function, Ei is the exponential-integral function, and the last line is valid in the 
vicinity of T^. For the class of models effectively described by the Morse potential, the fraction of bound pairs 
approaches zero linearly as T — > Tc, for any value of h (cf. Fig.^. The slope, as expected, reflects the particular 
choice of b. 



III. AN ALTERNATIVE VIEW OF THE TRANSITION: THERMAL STABILITY OF THE DOMAIN 

WALL 



A. Elementary dynamics 



The equation of motion corresponding to (2.1 



dV 

my„ = K{yn+i + 2/«-i - 2yn) - ^ — 

oyn 



(3.1) 



or, in the continuum limit, 



2/ = 



d'^y 1 dV 
dx^ m dy 



(3.2) 
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FIG. 5: The fraction of bound base pairs P{y < 2 A, T ) as a function of the T/Tc. 
soUd Une corresponds to the theoretical prediction (2.44). 



The triangles are TI values, whereas the 



where cq = cuoi. Relevant static configurations of the infinite chain with free ends, derivable from (3.2) are: 

(a) The uniform solution at the minimum of V{y), y{x) = 0; this configuration corresponds to the absolute minimum 
of the total pote ntial energy (equal to zero) and is stable with respect to small amplitude fluctuations; linearized 
solutions of ( ^.2[ ) are optical phonons. 

(b) For very large values of y, V{y) is almost a constant; therefore, expressions which verify (Py/dx^ — 0, i.e. of the 
form \iTaM^oc[y{x) — M] = C ■ x, where C > is an arbitrary constant, are a ppro ximate solutions in the sense that 



they correspond to local minima of the total potential energy; linearization of (3.2) around them leads to an acoustic 
phonon spectrum. 

(c) an exact, unbounded, domain-wall like solution 

1 



■In 



1 



(3.3) 



where d — i/^/2 R an d xo is an arbitrary constant. 

The solution ( |3.3[ ) is plotted in Fig. ^ for the case of the upper sign. It represents a nonlinear field configuration 
which " interpolates" from the stable minimum (a) to a particular member of the metastable configurations (b) , with a 
slope C = l/{ad), and hence equal contributions to the elastic and the on-site potential energy densities {D per site). 
In other words it is a profile of a double chain where the two strands stick together if x < xq, but strand separation 
grows linearly for all points x ^ xq; the energy of the solution contains a term which is proportional to the number 
of sites to the right of xq. More exactly, if lattice sites are numbered from to N, 



At zero temperature the profile (3_^) is not stable. Since 



E+^{xo±e) 



(3.4) 



(3.5) 



the wall will spontaneously move to the right, "zipping" back the unbound portion of the double chain. We will 
consider below how this instability changes under the infiucncc of temperature. 



B. Linearisation around the DW 



Consider small deviations with respect to (3.3), i.e 



X — xoje 



(3.6) 
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FIG. 6: The dashed hne represents the domain wall solution ( |3.3D as a func tion of the dimensionless space variable z. The full 
line represents the effective potential of the Schrodinger-like equation (3.7). The potential marginally fails to support bound 
states. Scattering states with eigenvalues less than the top of the potential (i.e. between dashed and dashed-double-dotted 
lines) are confined to the right-half of the line (low- frequency acoustic phonons); those with higher eigenvalues consist of both 
transmitted and reflected waves with different wavevectors and hence different dispersion relations in the two halves of the line 
(corresponding to both hard, i.e. optical, and soft, i.e. acoustic phonons, cf. text). 



where \oLj \ <^ a. The linearized eigenfunctions fj satisfy the equation 



1 ~ tanh z — sech z 



Rcu, 



(3.7) 







where z = {x — xq) / {2d) . Eq. (3/7) is Schrodinger-Uke with an effective potential drawn in Fig. |[ It can be shown jH 
that (3.7) has no bound states; scattering states are of two types: 



(I) if ujg < 2Rujq they are confined to the right half of the z-axis; for x ^ xq 



where 



fq cos [fc_z + Si] = cos [q{x — xq) + Si 



y = argF (ifc_) - argF - 1 + ^fc_ ) - argF 



2+ 



P = 2cj2/(i?cj§), k2 = 4 - fc2 , g = fc_/(2d) and therefore 



=ql ifqi< V2R 



(3.8) 
(3.9) 

(3.10) 



i.e. the "phonons" with the low frequencies are soft, acoustic phonons, as one expects them to be, since they correspond 
to atoms oscillating on the flat portion of the Morse potential. 

(II) in the high frequency case, i.e. if > 2Ruj^, scattering states extend over the whole line and have both 
transmitted and reflected components. Speciflcally, 



1/2 



where 



sinh ^ {k- — k+) 



sinhf (fc_ -I- fc+) 



(3.11) 



(3.12) 
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FIG. 7: Continuum limit dispersion relations for the soft (Eq. (|3.1C| ) as a dotted line) and hard (Eq. ( p.l3| ) as a solid line) 
phonons; they extend to ql — > oo. The discrete dispersion relations ( [4.l[ ) are also shown, with dashed line (hard) and dash- 
dotted line (soft). All curves have been drawn for R = 10.1; note the formation of a frequency gap between the two branches 
of the discrete spectrum for R > 2. 



T = 1 — p, fc^ = 4 — fc^, (/+ = k+/{2d) and therefore 



UJQ 



2R+iq+£y 



1/2 



if qe > V2R, a; > xo 
any q+, x xq 



(3.13) 



i.e. the dispersion relation is different in the soft (unbound) segment and the hard (bound) segment of the line (see 
Fig. 0). Formally, we can combine the two results for the dispersion into two branches, one optical and one acoustic, 
for each value of q, without any restrictions. However, it should be borne in mind that the acoustic phonons physically 
reside in the soft (unbound) and the optical ones in the hard (bound) segment of the line. 
For completeness we list the phase shifts 



SJj = 2argr 



2 + - (fc_ + k+) 



argF (1 + ik+) - argF (1 + ifc_) - tt 



S'jj = -2 argF 



2 + - (fc_ + fc^ 



2 argF 



2 + -{k^-k, 



2argF(l + ifc_) 



(3.14) 
(3.15) 



C. Free energy of the thermally dressed DW 



At finite temperatures, the domain wall is accompanied by a phonon cloud. The phonon cloud contributes three 
terms to the free energy, which arise from cases I (soft), II (soft), II (hard) presented in the previous subsection. 
Because these phonons reside in different - and in general unequal - portions of the chain, their contributions to the 
free energy depend on the position of the DW. Specifically, 

Fsoft = kBTY.^nphLj^'' (3.16) 

9 



Fhard = kBTY,^npnuj°P' (3.17) 

3 
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where contributions from soft modes have been combined. The sum over q can be replaced by a density of states 
determined by 

q{N-xo) = mr , n 0, 1, 2, ...(soft) (3.18) 
q{N + xo) = mr , n = 0, 1, 2, ...(hard) . (3.19) 

This results in a phonon-cloud free energy 

Fphonon cloud — ksTxQ / — In — — + ... (3.20) 

where the ellipsis denotes terms independent of xq and, since the integrand is ultraviolet-convergent, we have extended 
the upper limit of the integration to infinity. This is consistent with the continuum limit treated in this work; an 
ultraviolet divergence present in the ellipsis can be corrected by the introduction of a lattice cutoff but this is irrelevant 
for the purposes of the present argument. Introducing the correct dispersion relations for optical and acoustic phonons. 



we can evaluate [|14| the integral in (3.20). This results in a total free energy (DW plus phonon cloud) 




Fdw = const + ^/csT^^ - 2D j y . (3.21) 

Equation ( ^.2l| ) is the central result of this section. It describes in very simple terms why and when the phase 
transition occurs. At temperatures lower than 

2^/2D 2V2KD 
Tc = = 3.22 

the prefactor of xq is negative; the DW's natural tendency is towards high positive values of xq, i.e. it "zips" the 
system back to the bound configuration. Conversely, at temperatures higher than Tc, thermal stability is achieved by 
a high negative value of xq, i.e. the DW "opens up", unbinds the system. At the critical point, the DW is thermally 
(meta)stable at any position; this corresponds to physical realizations of the chain at critica lity. It should be noted 



that the value of Tc predicted by the above DW argument coincides with the exact result (2.14) of the continuum 
limit. This supports the validity of the alternative description of the phase transition in terms of the thermal stability 
of the DW. 



IV. DISCUSSION 



We have given a detailed account of the exact thermodynamics and the scattering function S{q, T) of a model used 
to describe thermal DNA denaturation, as well as other one-dimensional instabilities. The model deserves special 
attention for two reasons: (a) it is probably the simplest exactly solvable lattice model in one dimension which exhibits 
a true thermodynamic transition and satisfies all scaling laws (cf. below), and (b) the transition can be understood 
in terms of the thermal stability of a soliton-like nonlinear configuration. A number of comments are in order: 

1. None of the "prohibitions" of phase transitions in one dimension apply to this model. The theorem by Gursey [T^ 
and the van-Hove demonstration , extended by Ruelle ||l^ apply to systems with pair interactions only. The 
standard Landau argument, which covers systems with on-site potentials (e.g. 0*) rests on the finiteness of 
the DW energy. 



2. The critical exponents calculated analytically from the properties of the Schrodinger equation (2.10) satisfy 
all applicable scaling laws, i.e dv\\ = 2~a = '-f + 2(3 and 7 = (2 — Figs. 2-4 demonstrate very good 

agreement betw een the thermodynamic behavior calculated from numerical TI, and the one obtained exactly 



from Eq. ( |2lo| ). It appears that "universality" extends to prefactors and not just to exponents - provided that 



corrections due to the absolute value of critical point are taken into account. This in spite of the fact that the 
TI values refer to a system with a very high degree of discreteness {R = 10.1). 

3. It is instructive to view the phase diagram in the ih,T) plane (Fig. For positive values of the external field 
h, there is confinement at all temperatures. In the limit ft- ^ 0+, the system approaches a confined {T < Tc), 
or a deconfined (T > Tc) state. The TI (and the exact, Schrodinger) solutions presented here refer to the limit 
of the {h — 0~^,T) confined state as T — T~ (thick line). Finally, at any h < 0, the system is unstable. 
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FIG. 8: Free energy of the DW (in units o f D) for temperatures above (T = 'i/2Tc, solid line), below (T ; 
and at Tc (dotted line), according to Eq. (3.21); the irrelevant constant has been omitted. 



l/2Tc, dashed line) 



4. The scattering function (2.37) shows typical Ornstein-Zernike behavior. Its form demands that at least some 
of the weight cannot come from phonons. We will report details of critical dynamics in a separate paper. 
Preliminary results suggest the occurrence of a central peak. It will be interesting to relate critical dynamics 
with large-scale fluctuations of nonlinear configurations. 



5. The estimate ( 3.22 ) of the critical temperature can be extended to cover discrete systems, 
introducing the discrete dispersion relations 

2(1- cos((?^)) + 2R (optical) 
2(1 — cos{qi)) (acoustic) 



This is done by 



(4.1) 



m Eq. (|3!20|) which describes the free energy of the phonon cloud accompanying the DW. Numerical experiments 
have demonstrated |Q that the asymptotic slope of the DW is given by the same function of R for arbitrary 
levels of discreteness, i.e. lim„^oo(yn — yn-i)/i = l/{ad) = ^/2R/{a£) = \J2D j K j I. It follows that each site 
which finds itself in the high temperature phase contributes an elastic energy equal to (l/2)iir(?/„ — ?/„_i)^ = D 
and an on-site energy equal to D, i.e^ a total of 2D. This is a property of the DW, continuum or discrete. 



Therefore, the only modification to (3.21) and ( |3.22 ) comes from the entropy of the phonon cloud. In detail, 
using 



27r Jo 



dx ln[l — cos a; + JR] = In 



Vi? + \/i?T2' 



we obtain from ( 3.20 ) 



ksTxo 



phonon cloud 



[r(R)-r{o)] 



(4.2) 



(4.3) 



which results in the estimate 



Tr 



2D 



ks In 



yi?72 + Vl + i?/2 



(4.4) 



In the limit i? ^ 1, Eq. (4.4) reduces to Eq. ( 3.22| ). It can be seen in Fig. ( |lO| ) that the nonperturbative 
theoretical estimate (4.4) is in excellent agreement with TI numerical results even for values of R well into the 



14 



h 

▲ 



Confinement 

Instability 

► 




FIG. 9: Phase diagram in the (T, h) plane. At negative h, the system is unstable at all temperatures. 

discrete regime, R ^ 1. Systematic discrepancies appear at values of i? > 4 and seem to grow as R increases 
further. The origin of these discrepancies is currently under investigation; preliminary numerical results pc| ] 
suggest a significant complexity in the properties of the DW at values of i? > 4. 

Within the scope of the present work, we feel it is justified to state that (a) in the continuum approximation, and 
(b) even at moderate levels of discreteness, our understanding of a prototype one-dimensional phase transition 
can be considerably enhanced by making use of the concept of the thermal stability of a distinctly nonlinear 
entity. In more plain terms: formation of the DW can be thought of as "driving" the thermodynamic instability. 
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